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Abstract. The formal consideration of the concept of interac- 
tion in thermodynamic analysis makes it possible to deduce, in 
the broadest terms, new results related to the coevolution of inter- 
acting systems, irrespective of their distance from thermodynamic 
equilibrium. In this paper I prove the existence of privileged coevo- 
lution trajectories characterized by the minimum joint production 
of internal entropy a conclusion that extends Prigogine theorem 
to systems evolving far from thermodynamic equilibrium. Along 
these trajectories of minimum internal entropy production one of 
the system goes always ahead of the other with respect to cquilib- 



1. Introduction 

One of the primary objectives of non equilibrium thermodynamics is 
the analysis of open systems. As is well known, these systems maintain 
a continuous exchange (flow) of matter and energy with their environ- 
ment. Through these flows, open systems organize themselves in space 
and time. However, the behaviour of such systems differs considerably 
according to whether they are close or far from equilibrium. Close to 
equilibrium the phenomenological relationships which bind the flows to 
the conjugate forces responsible for them are roughly linear, i. e. of 
the type: 

Jt = ^LijXj, i, j = 1,2, . . .,n (1) 

3 

where Jj are the flows, Xj the generalized forces, and Ly are the so 
called phenomenological coefficients giving the Onsager Reciprocal Re- 
lations [6, 7]: 

L%j = Lji, % 7^ j (2) 

In these conditions, Prigogine's Theorem [8] asserts the existence of 
states characterized by minimum entropy production. Systems can 
assimilate their own fluctuations and be self-sustaining. 

l 
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Far from equilibrium, on the other hand, the phenomenology may 
become clearly non-linear, allowing the development of certain fluctu- 
ations that will reconfigure the system [1]. Thermodynamics of irre- 
versible processes can now provide information regarding the stability 
of the system but not on its evolution. Far from equilibrium no physi- 
cal potential exists capable of driving the evolution of systems [9] . The 
objective of the following discussion is just to derive certain formal 
conclusions related to that evolution in the case of interacting open 
systems evolving far from equilibrium. 



In the discussion that follow we will assume the following two basic 
assumptions: 

(1) There exist open systems whose available resources are limited. 

(2) Owing to these limitations, such systems compete with one an- 
other to maintain their necessary matter and energy flows. 

Both assumptions lead directly to the concept of interaction. But be- 
fore proposing a formal definition of that concept let us examine Figure 
1 in order to obtain an intuitive idea of the type of problem being dealt 
with. Figure 1 provides a schematic representation of three open sys- 
tems away from thermodynamic equilibrium. In the case (A) the sys- 
tem if not subjected to interaction with other systems, and maintains 
a through-flow of matter and energy which depends solely -without go- 
ing into phenomenological details- on its own degree of imbalance. In 
the case (B) the systems interact with each other, and their respective 
flows depend on the degree of imbalance of both systems simultane- 
ously. It is this situation which will be explored here as broadly as 
possible, particularly from the point of view of the coevolution of both 
systems. 

Let / and g be any two type C 2 functions (continuous functions 
with first and second derivatives which are also continuous functions) 
defined on the set K. of real numbers and such as: 



Function / will be referred to as the phenomenological function and g 
as the interaction function. The flow in a system will be given by J x = 
f(x) where x is the generalized force, a measure of the degree of the 
system's imbalance or distance from equilibrium. In consequence x > 
(x = at equilibrium); x will also be used to designate the system. 
Let us now justify the above constraints (3)- (5) on / and g. Firstly, 
/(0) = g(0) = indicates that in the absence of generalized forces 
(thermodynamic equilibrium) there is neither flows nor interactions. 



2. Definitions 



/ and g are strictly increasing, and /(0) = g(0) = 
(/ — g) is strictly increasing 
(/ — g)' is increasing. 



(3) 
(4) 
(5) 
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Sink 



Figure 1: Circles represent open systems far from thermodynamic equilibrium. 
Rectangles represent significant parts of their environment (resources and sinks). 
Matter and energy flows are shown as large arrows crossing systems from resources 
to sinks. Yellow vertical arrows x and y represent the generalized forces responsible 
for the Bows. Curved arrows g(x) and g(y) represent the interactions between the 
systems. In the absence of interaction with other systems, Bows depend only on 
the intensity of the own force that maintain the system away from equilibrium 
(Jx = f( x ))- When two systems compete with one another in order to maintain 
their Bows, these Bows depend on the intensity of the generalized forces (degree 
of imbalance) in both systems simultaneously. In these conditions, the Bows are 
expressed by expression of the type J x = f(x) — g(y); J y = f(y) — g(x) (see text). 



The increasing nature of both functions indicate that the intensity of 
flows and interactions increase with the generalized forces, or in other 
words, as we move away from equilibrium. Constraints (4) and (5) 
imply that systems are progressively as sensitive to their own forces as 
they are to the forces of the other interacting system, or more so. 

Given two systems with the same phenomenological function / and 
forces x and y respectively, an interaction can be said to exist between 
them if their respective flows can be described by the following expres- 
sions: 

J x = f(x) + C xy g(y), -1 < C xy < 1 (6) 
J y = f{y) + C yx g{x), -1 < C yx < 1 (7) 

where C xy and C yx are the interaction coefficients. We shall examine 
here the double negative, the positive-negative, and the double positive 
interaction, assuming -while still speaking in general terms- that C xy = 
C yx = ±1. Or in other words, that: 

(j x = f(x)±g(y) 

1 Jy = f(y) ± g(x) 



In addition, the only (x, y) pairs permitted will be those which give a 
positive value for expressions (8) above; (x, y) pairs giving f(x)—g(y) = 
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shall be termed the points of extinction of system x. The same applies 
to system y 

3. Entropy production 

As is well known, the entropy balance for an open system can be ex- 
pressed as: 

dS = diS + d e S (9) 

where d^S represents the entropy production within the system due to 
flows, and d e S is the entropy exchanged with its surroundings. The 
second law of thermodynamics dictates that always d^S > (zero only 
at equilibrium). 

One of the most interesting aspects of Thermodynamics of Irre- 
versible Processes is the inclusion of time in its equations: 

^to=Si = Yl x j J i ( 10 ) 

3 

In our case, for system x we have: 

S ix =x[f(x)-g(y)]>0 (11) 



dt 

and for system y: 

din S 



^ r = S iy = y[f(y)-g(x)]>0 (12) 



and the joint production of internal entropy: 

di:r-S disr.S 



S(x, y) = x[f(x) - g{y)\ + y[f(y) - g(x)} (13) 

Graphically S(x, y) is a surface in the space defined by the axes S, 
X and Y (Figure 2). A curve in the plane XY represents a possible 
coevolution history for both systems (as long as x and y satisfy (11) 
and (12)). Consequently, this type of curves will be termed as coevolu- 
tion trajectories. The projection on the surface S(x, y) of a coevolution 
trajectory represents its cost in terms of internal entropy production. 
As we will immediately see, there exist special coevolution trajectories 
characterized by the minimum entropy production inside the systems 
evolving along them. They will be referred to as trajectories of mini- 
mum entropy (TME). Since the minimal entropy production inside the 
systems means the maximum regularity in their spacetime configura- 
tions, the points of a TME represents, states of maximum stability in 
systems interacting far away from thermodynamic equilibrium. 

We will now analyze surface S(x, y) for negative-negative, positive- 
negative and positive-positiveinteractions. And we will do it following 
a common mathematical method which consists in examining the in- 
tersections of the surface with one or more families of planes (Figure 
2). In our case, three families will be used: the family of planes parallel 
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Figure 2: Left: Curve 1-2 represents any coevolutionary history between systems 
x and y. Curve l'-2', its projection on S(x,y), represents the joint production of 
internal entropy throughout that history. Right: Intersections of surface S(x,y) 
with the family of planes parallel to SX. 



to the plane SX; the family of planes parallel to the plane SY; and 
the family of planes of the form y = c — x, which are perpendicular to 
the bisector y = x. 

3.1. Negative-negative interaction. According to (8), in the case 
of a negative- negative interaction we will have: 



J X 



/0) -g(y) 
f(y) - a(x) 



(14) 



In these conditions, surface S(x, y) is given by: 

S(x, y) = x[f(x) - g(y)} + y[f(y) - g(x)} (15) 
Let us analyze its intersections with our three family of planes. 



Planes parallel to SX 

Let b be any real positive number. Consider the plane y = b. Its 
intersection with S(x, y) will be: 

S(x, b) = h(x) b = x[f(x) - g(b)} + b[f(b) - g(x)} (16) 

If we derive h(x)b we will have: 

= h\x) b = f(x) - g(b) + xf(x) - bg'(x) (17) 

According to restrictions (3), (4) and (5), h'(x)b is strictly increasing, 
and according to the same restrictions: 

^(0)6 = -^(6) - 6^(0) < (18) 

^(6) 6 = f(b) - g{b) + b[f'(b) - g\b)\ > (19) 
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Figure 3: Intersection of surface S(x,y) with the plane y = b parallel to the plane 
SX. 

Therefore, and in accordance with Bolzano's theorem, there is a point 
Xb such that: 

h'{x b ) b = 0; < x b < b (20) 
And, since h'{x)b is strictly increasing, we have: 

ti\x) b > 0; Vx > (21) 
Consequently Xb is a minimum of S(x, b) (Figure 3) 

So, the intersection of surface S(x,y) with the plane y = b parallel to 
the plane SX is a curve h(x)b with a minimum. Or in other terms, for 
every y > there is a value of x that minimizes the joint production 
of internal entropy in both systems. The couples (x, y) satisfying this 
condition define a curve T\ (x, y) in the plane XY whose equation is: 

f(x) - g(y) + xf(x) - yg'(x) = (22) 

which results from the fact that h\x) y = for each (x, y) in Ti(x,y). 
Consequently, the couples (x, y) satisfying (22) form a coevolution tra- 
jectory characterized by the minimal joint production of entropy inside 
the systems evolving along its points, that is to say a TME. According 
to (20), it holds the following St. Matthew inequality: 

V(x,y) E Ti(x, y) : x < y (23) 

Note that according to the above reasoning this asymmetry is universal, 
independent of the particular functions / and g, provided they verify 
constraints (3)-(5). 

Planes parallel to SY 

The same reasoning above now applied to the intersections of S(x, y) 
with planes parallel to SY leads to a new coevolution trajectory T 2 (x, y) 
in the plane XY whose equation is: 

f(y)-g(x)+yf'(y)-xg'(y)=0 (24) 
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Figure 4: Symmetry of the trajectories of minimum entropy production. 

and whose points are also characterized by the minimal joint production 
of internal entropy in the systems evolving along them, i.e. A new 
TME. As in the case of T\(x,y), and for the same reasons, a new St. 
Matthew inequality holds: 

\f(x,y) E T 2 (x,y) : y < x (25) 

It is easy to see that Ti(x,y) and T 2 (x,y) are symmetrical with 
respect to the bisector y = x. In fact, for y — equation (22) becomes: 

f(x) + xf(x) = (26) 

which only holds at x — 0. The point (0, 0) belongs, therefore, to 
the trajectory. The same applies to (24). Consequently, the point 
(0,0), which corresponds to thermodynamic equilibrium, belongs to 
both trajectories. On the other hand, it is evident that: 

T 1 (x,y) = T 2 (y,x) (27) 

Therefore, both trajectories are symmetrical with respect to the bisec- 
trix y = x. 

Intersections with planes of the form y = c — x 

Let c be any positive real number, the intersection of S(x, y) with the 
plane y = c — x will be: 

h(x) c _ x = x[f(x) - g(c - x)\ + (c - x) [f(c - x) - g(x)] (28) 

whose derivative is 

h{x)' c _ x = f(x) - g(c -x) + xf'(x) + xg'(c - x) 

- f(c -x)+ g{x) - (c - x)f(c -x)-(c- x)g'(x) (29) 
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This derivative vanishes at point x = c/2: 




= (30) 



And according to restrictions (3), (4): 



h "Q= i (f'® + s'(d) + < f "®s"(ii)) >0 m 

we conclude that point (c/2, c/2) is a minimum of the intersection 
h(x) c _ x . The same can be said of each point in the bisector y = x. In 
consequence, this line is a new trajectory of minimum entropy. It will 
be referred to as T 3 (x,y). 

s 
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Figure 5: S(x, y) plotted in four sections orthogonal two by two. Each of those 
sections results from the intersection of S(x, y) with the planes x = 0; y = 0; x = a, 
y = a. A fifth section resulting from the intersection of S(x, y) with the plane 
y = a — x has also been plotted. 

The above results allow us to depict surface S(x, y) in the space defined 
by the axis X, Y and S (Figure 5). Although the bisector y = x is also 
a TME, it is not the most efficient in terms of entropy production. In 
effect, for each point (6, b) of this line, two points (x, b) and (b, y) exist 
at which the interacting systems produce the less possible amount of 
internal entropy. The first of these points belong to Ti(x, y), the second 
to T 2 (x,y) (Figure 6). 

As we have just seen, for each point (a, b) of a coevolution trajectory, 
including the points of T 3 (x,y), there is a point (x, b) in T±(x,y) and 
a point (a,y) in T 2 (x,y) such that the joint entropy production inside 
the systems is the less possible one. This evidently applies to any 
point (c, d) in the plane XY. In addition, and as a consequence of St. 
Matthew inequalities, one of the systems is always ahead of the other 
with respect to equilibrium. It holds therefore the following: 
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Figure 6: Although the bisector y — x is also a TME, it is not the most efficient 
one in terms of internal entropy production. 

Theorem 1. For every negative-negative interaction two trajectories 
of minimum entropy exist formed by a succession of states character- 
ized by the minimum entropy production within the interacting systems. 
Furthermore, one of the systems is always ahead of the other with re- 
spect to thermodynamic equilibrium. 



3.2. Positive-negative interaction. According to equations (6)-(7), 
and under the same restrictions of the negative-negative case, we will 
have now: 



The same considerations on the joint production of internal entropy 
we made in the case of the double negative interaction allow us now to 
express the entropy production in the positive-negative case as: 



production inside the interacting systems. We will analyze it with the 
same intersection method. 



J x = f(x)+g(y) 

J v = f{y)-g{x) 



(32) 




Intersections parallel to SX 




(34) 
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And its derivative: 

h'(x) b = f(x) + g{b) + xf'(x) - bg'(x) (35) 

which is an strictly increasing function in accordance with the restric- 
tions (3), (4) and (5). In addition, for x = equation (35) becomes: 

h'(0) b = g(b) - bg'(0) (36) 

While for x = b: 

ti{b) b = fib) + g(b) + bifib) - g\b)) > (37) 

Consequently, and according again to Bolzano's theorem, if: 

g(b) - bg'(0) < (38) 

there will be a point x b such that: 

h'(x b ) b = (39) 

Taking into account that h'(x) b is strictly increasing, we have: 

h"(x) b > 0,Vx > (40) 

Therefore x b will be a minimum. Let us define gi(x) as g(x) — xg'{0). 
According to (38) there will be a trajectory of minimum entropy T±(x, y) 
if: 

gt(x) < 0, Vx > (41) 

As in the case of the negative-negative interaction, its equation, will 
be: 

f[x) + g{y) + xf'(x) - yg'ix) = (42) 



Intersections parallel to SY 

Let a be any positive real number. The intersection of (33) with the 
plane x = a is: 

h(y) a = aifia) + g{y)) + y(f(y) - g{a)) (43) 
and its derivative: 

h\y) a = ag'(y) + f{y) - g{a) + yf'(y) (44) 

which is a strictly increasing function according to restrictions (3), (4) 
and (5). For y = equation (44) becomes: 

h'i0) a = ag'i0)-gia) (45) 

while for y = a we have: 

h'(a) a = ifia) - gia)) + a(g'(a) + /'(a)) > (46) 

Therefore if: 

ag'{0)-g{a)<0 (47) 
there will be a point y a such that: 

h'iy a ) a = (48) 
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Taking into account the strictly increasing nature of h'(y) a we will have: 

h"(y) a > 0, Vy > (49) 

Therefore y a will be a minimum. Let us define g 2 (x) as xg'(0) — 
g{x). According to (47) there will be a trajectory of minimum entropy 
Ti(x,y) if: 

g 2 {x) < 0, Vx > (50) 

Its equation will be: 

f(y)-g(x)+yf'(y)+xg'(y) = (51) 

Since 

gi (x)=g(x)-xg'(0) (52) 
g 2 (x)=xg'(0)-g(x) (53) 

it is clear that: 

0i (a) = ~92(x) (54) 

Both functions will only coincide at point (0, 0). In addition, if g±(x) < 
then g 2 (x) > 0, and viceversa. Therefore, if g±(x) < there exists 
Ti(x,y), but not T 2 (x,y). Similarly, if g 2 (x) < then there exists 
T 2 (x,y), but not Tx(x,y). And since xg'(0) and g{x) are two strictly 
increasing functions growing from (0,0), it must exclusively hold one 
of those two alternatives; and then only one trajectory or minimum 
entropy, either Ti(x,y) or T 2 (x,y), will exist. 



s 




i 

X=Y 



Figure 7: In the case of the positive-negative interaction, only a trajectory of 
minimum entropy exists. 



Intersections with planes of the form y = c — x 
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Being c any positive real number, the intersection of S(x, y) with the 
plane y = c — x will be: 

h(x) c _ x = x[f(x) + g(c - x)] + (c-x) [f(c - x)) - g{x)\ (55) 

and its derivative: 

h'(x) c _ x = f(x) + g(c - x) + xf(x) - xg'(c - x) 

— f(c — x) + g(x) — (c — x)f'(c — x) — (c — x)g\x) (56) 



For x = c/2 we have: 



/( 



c 
2 



2 (^ C ^ 2° ' ^ C 
~ 51 V2/ 2 9 \2, 

Accordingly, the point (c/2, c/2), and any other of the bisector y = x, 
will be a maximum or a minimum if: 

5-(x) = xg'(x) (57) 

Point (c/2, c/2) will be a minimum if: 

^(|) = ^(f) + c(r(|) > (58) 

In this case the bisector ?/ = i is a TME. It will be referred to as 
T 3 (x, y). We have, then, three alternatives for the existence of TME in 
the case of positive-negative interactions: 

Ti(x,y) exists if: xg'(0) > g(x) (59) 

T 2 (x,y) exists if: xg'(O) < g(x) (60) 

T 3 (x,y) exists if: xg'(x) = g(x) (61) 

T 3 (x,y) also requires condition (58). 

As in the case of the double negative interaction, given a point (a, b) of a 
coevolutionary trajectory, there exist either a point (x, b) in Ti(x, y) or 
a point (a, ?/) in T 2 (x, y) such that the interacting systems produce the 
less possible amount of entropy inside the systems. Unlike the double 
negative interaction, in the positive-negative case only one TME exists. 
If that trajectory is: 

f(x) + g(y) + xf'(x) - yg'(x) = (62) 

it holds: 

y f'(x) f'(x) f'(x) [b6) 
According to restrictions (3)- (4)- (5), the first fraction is positive and 
less than 1, while the second and the third ones are positives. Therefore 
x < y, and system y is always ahead of system x. 
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The other possible TME is: 

f(y)-g(x)+yf'(y)+xg'(y) = (64) 

whence: 

_ g(x) - f(y) - xg'jy) 

y ~ m ( J 

In accordance with restrictions (3)-(4)-(5), if y > x the numerator 
would be negative, and being the denominator always positive, we 
would have a negative value for y, which is impossible (forces x and 
y are always positive). So y is less than x. Consequently system x is 
always ahead of system y. It holds, therefore, the following: 

Theorem 2. For every positive-negative interaction a trajectory of 
minimum entropy exists formed by a succession of states character- 
ized by the minimum entropy production within the interacting systems. 
Furthermore, one of the systems is always ahead of the other with re- 
spect to thermodynamic equilibrium. 

3.3. Positive-positive interaction. In the case of a double positive 
interaction, in the place of (8) we will have: 



J x = f{x) + g{y) 
J y = f(y) + g(x) 



(66) 



and instead of (15): 

S(x, y) = x[f(x) + g{y)} + y[f{y) + g{x)} (67) 



Whence: 



= f(x) + g(b) + Xf\x) + bg'(x) (61 

dS(a,y) 



dy 



f(y) + g(a) + Yf'(y) + ag'(y) (69) 



Both derivatives are positive, strictly increasing and only vanish at 
point (0, 0) corresponding to thermodynamic equilibrium. Intersections 
S(x,b), and for the same reasons intersections S(a,y), have neither 
maximums nor minimums. 

The intersection of surface S(x, y) with the plane y = c — x is now:: 

h(x) c _ x = x(f(x) + g(c - x)) + (c - x)(f(c - x) + g(x) (70) 
and its derivative: 



h'(x) c - x = f(x) + g(c -x)+ xf'(x) - xg'(c - x) 

- f(c -x) - g(x) - (c - x)f(c -x) + (c- x)g'(x) (71) 
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which vanishes at x = c/2): 



2/ c -x \2/ \2/ 2" V2/ 2" V2 

Therefore, there is a minimum at x — c/2 if the second derivative 
h"(x) c - x is positive. That is to say if: 

2(f'(x) - g'(x)) + x(f"(x) + g"(x)) > (72) 

In short, the only possible TME in the case of a double positive in- 
teraction is the bisectrix y = x. This conclusion together with the 
above two St. Matthew asymmetries allow us to state that stability in 
interacting systems requires either asymmetry or equality depending 
on the interaction nature: asymmetry for competition and equality for 
cooperation. 



4. St. Matthew Theorem 

'St Matthew Effect' is the title of a paper by Robert k. Merton pub- 
lished in Science in the year 1968. The main objective of that paper 
was the analysis of the reward and communication systems of science. 
In its second page we can read: "... eminent scientists get dispropor- 
tionately great credit for their contributions to science while relatively 
unknown scientists tend to get disproportionately little credit for com- 
parable contribution". [5]. Or in other more general terms: the more 
you have the more you will be given, a pragmatic version of the so 
called St Matthew principle: "For unto every one that hath shall be 
given, and he shall have abundance, but from him that hath not shall 
be taken away even that which he hath". (St. Matthew 25:29). So- 
ciologists, economists, ethologists and evolutionary biologists, among 
many others, had the opportunity to confirm the persistence of this St. 
Matthew asymmetry, which invariably appears in every conflict where 
any type of resource, including information, is at stake. 

We have just proved, in the broadest terms, the existence of asym- 
metric coevolution trajectories for systems that interact according to 
certain formal definitions. These trajectories are thermodynamically 
relevant because their points represent states very far from equilibrium 
characterized by the minimum joint production of internal entropy in 
the interacting systems. Along these trajectories one of the system is 
always ahead of the other with respect to equilibrium. It has also been 
proved that those coevolution trajectories of minimum entropy exist if, 
and only if, at least one of the interactions is negative in the formal 
sense defined above. It has therefore been proved the following: 



St Matthew theorem. For every binary interaction in which at least 




one of the interactions is negative, a coevolution trajectory exists char- 
acterized by the minimum joint production of internal entropy within 
the systems evolving along it, and such that one of the systems always 
goes ahead of the other with respect to thermodynamic equilibrium. 

St Matthew theorem states, therefore, the same conclusion for systems 
far from thermodynamic equilibrium as Prigogine theorem for systems 
close to it. 



5. Generalization 

If we consider n systems instead of 2, the following expressions would 
be obtained for flows: 

J Xi = f(xi) + ^2C ij g(x j ), -1 < Ctj < 1, i,j = 1,2, . . .n (73) 

3 

where C^- represents (as the binary case) the type and degree of in- 
teraction between systems X{ and Xj. The joint production of internal 
entropy could thus be expressed as: 

(74) 

where S(x±,X2, ■ ■ ■ x n ) is a hypersurface on M n+1 in which we can con- 
sider ^2{n — 1) three-dimensional subspaces SxiXj in order to analyze 
the interaction between the system and the system xj using the same 
method as in the binary case. 



S(xi,x 2 , 



f( x i) + ^2 C ij9(xj] 
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6. Discussion 

The most relevant feature of an open system is its ability to exchange 
matter and energy with its environment. But not all processes involved 
in those exchanges are equally efficient. Far form equilibrium, systems 
are subjected to the non linear dynamic of fluctuations. In those con- 
ditions, the most efficient process are those that generate the lower 
level of internal entropy because the excess of entropy may promote 
an excess of fluctuations driving the system towards instability. We 
have proved the existence of coevolution trajectories whose most re- 
markable characteristic is just the minimum level of internal entropy 
production. The states of the systems evolving along those trajectories 
are, therefore, the most efficient ones in terms of self maintenance. In 
these conditions and taking into account the tremendous competence 
and selective pressure suffered by many systems, as is the case of the 
biological or the economical systems, the trajectories of minimum en- 
tropy should be taken as significant references. They in fact represent 
histories of maximum stability in open systems interacting very far 
from thermodynamic equilibrium, which is very a common situation in 
the real world. Apart from its own existence, it is remarkable the asym- 
metrical way the systems evolve along them (St Matthew theorem). A 
result compatible with the empirical observations of evolutionary biol- 
ogy that biologists know long time ago and that is usually referred to as 
St Matthew principle [4] . It is also confirmed by a recent experimental 
research related to recursive interactions by means of logistic functions 
[3]. 
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